function yyy=equation_difference(xxx,params)

gammaI=xxx(1);
xI=xxx(2);
A=xxx(3);
gammaE=xxx(4);
xE=xxx(5);
B=xxx(6);
gq=xxx(7);
xe=xxx(8);
gammae=xxx(9);
C=xxx(10);

alpha=params.alpha;
pi=params.pi;
eta=params.eta;
r=params.r;
lambdaE=params.lambdaE;
lambdaS=params.lambdaS;
lambdaI=params.lambdaI;
lambdaC=params.lambdaC;
lambdae=params.lambdae;
f=params.f;
xiE=params.xiE;
xiI=params.xiI;
xie=params.xie;

% equation (1)
gammaI_est=xiI*(alpha/(1-alpha))*xI^(1/(1-alpha))*pi;
eq1=abs(gammaI_est-gammaI);

% equation (2)
TEMPTEMP=((1-lambdaC)^(eta-1)+(1+lambdaI)^(eta-1)-1);
TEMP=(((1-alpha)/xiI)*(A/pi*TEMPTEMP+f*B/pi));
xI_est=(TEMP)^((1-alpha)/alpha);
eq2=abs(xI_est-xI);

tau=xE+xe;
psi=xI+tau;

% equation (3)
TEMP=(1-lambdaS)^(eta-1)-1;
A_est=(pi+gammaI)/(r-(2-eta)*gq-tau*TEMP+psi);
eq3=abs(A_est-A);

% equation (4)
gammaE_est=xiE*(alpha/(1-alpha))*xE^(1/(1-alpha))*pi;
eq4=abs(gammaE_est-gammaE);

% equation (5)
TEMP=((1-alpha)/xiE)*(A*(1+lambdaE)^(eta-1)/pi+B/pi);
xE_est=(TEMP)^((1-alpha)/alpha);
eq5=abs(xE_est-xE);

% equation (6)
B_est=gammaE/(r-(2-eta)*gq+psi);
eq6=abs(B_est-B);

% equation (7)
TEMP1=((1+lambdaE)^(eta-1)+(1-lambdaS)^(eta-1)-2)/(eta-1);
TEMP2=((1+lambdaI)^(eta-1)+(1-lambdaC)^(eta-1)-2)/(eta-1);
TEMP3=((1+lambdae)^(eta-1)+(1-lambdaS)^(eta-1)-2)/(eta-1);
gq_est=xE*TEMP1+xI*TEMP2+xe*TEMP3;
eq7=abs(gq_est-gq);

% equation (8)
TEMP=A*(1+lambdae)^(eta-1)/pi+B/pi-C/pi;
xe_est=((1-alpha)/xie*TEMP)^((1-alpha)/alpha);
eq8=abs(xe_est-xe);

% equation (9)
gammae_est=xie*alpha/(1-alpha)*(xe)^(1/(1-alpha))*pi;
eq9=abs(gammae_est-gammae);

% equation (10)
C_est=gammae/(r-(2-eta)*gq);
eq10=abs(C_est-C);

eq = [eq1^2 eq2^2 eq3^2 eq4^2 eq5^2 eq6^2 eq7^2 eq8^2 eq9^2 eq10^2];

yyy = sum(eq);
end

